Search for: All records

Abstract In this article, we study the spectral volume (SV) methods for scalar hyperbolic conservation laws with a class of subdivision points under the Petrov–Galerkin framework. Due to the strong connection between the DG method and the SV method with the appropriate choice of the subdivision points, it is natural to analyze the SV method in the Galerkin form and derive the analogous theoretical results as in the DG method. This article considers a class of SV methods, whose subdivision points are the zeros of a specific polynomial with a parameter in it. Properties of the piecewise constant functions under this subdivision, including the orthogonality between the trial solution space and test function space, are provided. With the aid of these properties, we are able to derive the energy stability, optimal a priori error estimates of SV methods with arbitrary high order accuracy. We also study the superconvergence of the numerical solution with the correction function technique, and show the order of superconvergence would be different with different choices of the subdivision points. In the numerical experiments, by choosing different parameters in the SV method, the theoretical findings are confirmed by the numerical results. 

In this paper, we consider a class of discontinuous Galerkin (DG) methods for one-dimensional nonlocal diffusion (ND) problems. The nonlocal models, which are integral equations, are widely used in describing many physical phenomena with long-range interactions. The ND problem is the nonlocal analog of the classic diffusion problem, and as the interaction radius (horizon) vanishes, then the nonlocality disappears and the ND problem converges to the classic diffusion problem. Under certain conditions, the exact solution to the ND problem may exhibit discontinuities, setting it apart from the classic diffusion problem. Since the DG method shows its great advantages in resolving problems with discontinuities in computational fluid dynamics over the past several decades, it is natural to adopt the DG method to compute the ND problems. Based on [Q. Du, L. Ju, J. Lu and X. Tian,Commun. Appl. Math. Comput. 2 (2020) 31–55], we develop the DG methods with different penalty terms, ensuring that the proposed DG methods have local counterparts as the horizon vanishes. This indicates the proposed methods will converge to the existing DG schemes as the horizon vanishes, which is crucial for achievingasymptotic compatibility. Rigorous proofs are provided to demonstrate the stability, error estimates, and asymptotic compatibility of the proposed DG schemes. To observe the effect of the nonlocal diffusion, we also consider the time-dependent convection–diffusion problems with nonlocal diffusion. We conduct several numerical experiments, including accuracy tests and Burgers’ equation with nonlocal diffusion, and various horizons are taken to show the good performance of the proposed algorithm and validate the theoretical findings. 

In this paper, we study the central discontinuous Galerkin (DG) method on overlapping meshes for second order wave equations. We consider the first order hyperbolic system, which is equivalent to the second order scalar equation, and construct the corresponding central DG scheme. We then provide the stability analysis and the optimal error estimates for the proposed central DG scheme for one- and multi-dimensional cases with piecewise P k elements. The optimal error estimates are valid for uniform Cartesian meshes and polynomials of arbitrary degree k  ≥ 0. In particular, we adopt the techniques in Liu et al . ( SIAM J. Numer. Anal. 56 (2018) 520–541; ESAIM: M2AN 54 (2020) 705–726) and obtain the local projection that is crucial in deriving the optimal order of convergence. The construction of the projection here is more challenging since the unknowns are highly coupled in the proposed scheme. Dispersion analysis is performed on the proposed scheme for one dimensional problems, indicating that the numerical solution with P 1 elements reaches its minimum with a suitable parameter in the dissipation term. Several numerical examples including accuracy tests and long time simulation are presented to validate the theoretical results.